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Abstract - Recently, a new quantum model - two-dimensional generalization of the Scarf II - 
was completely solved analytically by SUSY method for the integer values of parameter. Now, 
the same integrable model, but with arbitrary values of parameter, will be studied by means 
of supersymmetrical intertwining relations. The Hamiltonian does not allow the conventional 
separation of variables, but the supercharge operator does allow, leading to the partial solvability 
of the model. This approach, which can be called as the first variant of SUSY-separation, together 
with shape invariance of the model, provides analytical calculation of the part of spectrum and 
corresponding wave functions (quasi-exact-solvability). The model is shown to obey two different 
variants of shape invariance which can be combined effectively in construction of energy levels and 
wave functions. 



Introduction. — Supersymmetrical Quantum Me- 
chanics gave a new impetus to development of new an- 
alytical methods in study of different quantum models [1] 
- [3] . The main ingredients and results of SUSY approach 
in Quantum Mechanics are isospectral (or, almost isospec- 
tral) quantum systems and supersymmetrical intertwining 
relations [1] , shape invariance [I] - [5] , higher order super- 
charges [5] , supersymmetry of multidimensional [TU] - [Tl] , 
multiparticle [15) and matrix |16j quantum systems. 

One of the directions of research within the framework 
of this approach is investigation of two-dimensional quan- 
tum models which are not amenable to standard |17j sep- 
aration of variables. Recently, two different variants of su- 
persymmetrical (SUSY) separation of variables were pro- 
posed [8] , [18] - [20] for analysis of such two-dimensional 
quantum systems. Both procedures are based on super- 
symmetrical intertwining relations [T] , |10] with second or- 
der supercharges [13] - [14], [18] and shape invariance [4], 
0, [8], [18], (BJ, [7] property. The first variant is applica- 
ble if the intertwining operator (supercharge) allows the 
standard separation of variables as it happens for super- 
charges with Lorentz metric. It was used in practice for 
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two-dimensional generalizations of Morse [5] , [T5] , Poschl- 
Tellcr [21], and periodic Lame [TH] models. The second 
variant works if one of the intertwined Hamiltonians allows 
the standard separation due to specific choice of param- 
eters, but its supcrpartner still does not. Such situation 
gives a chance to solve the problem completely - to find en- 
ergies and wave functions of all bound states. This proce- 
dure was applied recently to Morse [5D] , Poschl- Teller [33] 
models, and quite recently [23] , to the new model - two- 
dimensional generalization of Scarf II model. In the latter 
case, one of intertwined Hamiltonians is separable if one 
of parameters takes only negative integer values a = —k. 

In the present paper, the same two-dimensional Scarf 
II model but with arbitrary values of parameter a will be 
studied by means of the first procedure of SUSY separa- 
tion of variables. This method provides a part of energy 
spectrum and corresponding wave functions in analytical 
form. Thus, the model belongs to the class of partially (i.e. 
quasi-cxactly-solvable) models - the intermediate class be- 
tween completely (exactly solvable) and analytically un- 
solvablc models. This class was considered starting from 
80-ties [24] - [28]. In particular, the elegant algebraic 
method of construction of one-dimensional quasi-exactly- 
solvable (and sometimes, of exactly solvable) quantum 
models was elaborated in [25] - [27] . This approach is ap- 
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plicable to two-dimensional problems as well, but only in 
curved spaces with nontrivial metrics |26j . The approach 
of the present paper allows to study the different class of 
quantum models which do not allow standard separation 
of variables. It is necessary to notice that by construction 
these models are intcgrable with the symmetry operators 
of fourth order in momenta. In Section 2, the zero modes 
of supercharge will be found, and their linear combina- 
tions, which are the eigenfunctions of Hamiltonian, will 
be built. In Section 3, the shape invariance of the model 
will be used to enlarge the variety of known wave func- 
tions. It will be shown that the second shape invariance 
exists for this model leading to additional relations be- 
tween different wave functions. 

Wave functions in subspace of zero modes of Q + . 

Formulation of the model. The main tools of the 
supersymmetrical approach in two-dimensional Quantum 
Mechanics are the supersymmetrical intertwining rela- 
tions: 



ff (1) Q + = Q+ff (2) ; Q-HW=HWQ- 



(1) 



for two partner two-dimensional Hamiltonians of 
Schrbdinger form 



H (i) = _AP9+yW(2); i 



1,2; 



x = (x 1 ,x 2 )(2) 



_d_ 
dxi 



with mutually conjugate supercharges Q of second or- 
der in derivatives. First of all, one has to find solutions of 
Eq. ((TJ, i.e. to find such potentials V^ 1,2 \x) and such coef- 
ficient functions of second order supercharges Q ± that (JTJ 
are fulfilled. In terms of these unknown functions, we deal 
with a complicate system of nonlinear differential equa- 
tions of second order. Due to suitable choice of ansatzes, 
a list of particular solutions of (JTJ was found [13] . |18j . 
A part of these solutions was shown to allow the analyti- 
cal construction of spectra and wave functions: depending 
on chosen values of parameters, the partial [5], [HJ, [T5] . 
[TB] and/or complete [5D] , [22] solutions of corresponding 
models were obtained. 

Recently, new results [6] obtained in one-dimensional 
shape invariance allowed to find also new solutions [7] of 
{!]) in two-dimensional framework. In particular, the two- 
dimensional generalization of Scarf II model was among 
these new solutions, and just this system will be considered 
below. The potentials 

y(i),(2) = _2A 2 a(a T l)( i . - i . ) + 

i 2B(A+\)sinh(2\x 1 ) + (B 2 -A 2 -2A\) . 



2B{A+X) sinh{2Xx 2 ) + (B 2 -A 2 -2AX) 
cosh 2 (2Aa:2) ' 



(3) 



(4) 



and the second order supercharges: 

Q+ = (Q-)t = Ad+d- + 4aAtanh(Ax+)d_ + 
+4aAcoth(Ax_)<9+ + 4a 2 A 2 tanh(Ax+) coth(Ax_) 

2B(A+X) sinh(2Xx 1 ) + (B 2 -A 2 -2AX) 

cosh 2 (2Aici) '~ 

2B(A+X)sinh(2Xx 2 ) + {B 2 -A 2 -2AX) 
' cosh 7 (2AaJ7) ' 

solve the intertwining relations (pQ)(x± = x\ ±X2,<9± = 
d/dx±; X, a are real parameters, and A, B > 0). As is typ- 
ical to the approach, both potentials ([3]) are not amenable 
to standard separation of variables, but they correspond 
to integrable Hamiltonians ((2|) with symmetry operators of 
fourth order in derivatives: R^ — Q + Q~ , R^ = Q~Q + . 
Quite recently, it was proven that this model is com- 
pletely (exactly) solvable for the values of parameter a = 
—k. The whole spectrum of bound states and wave func- 
tions were found analytically by means of intertwining re- 
lations (TTJ) and shape invariance of ([3]) under the change 
a — > a — 1 (see Section 3 below) . 

Zero modes of supercharge Q + . The same model will 
be studied below for arbitrary values of a in the framework 
of SUSY separation of variables of the first kind (see [8] , 
|21j . }19j . (18)). The basic idea is that the supercharge 
Q + for Lorentz (hyperbolic) form of metric, like in (01, is 
amenable to conventional separation of variables [5] , [TS] . 
Indeed, after a suitable similarity transformation: 

Q+ = exp (x)q + exp (-%); 
x{x) = — a In | cosh(Ax+) sinh(Ax_)|, 
we obtain the operator with separated variables: 
q+ = dl-d 2 2 -f{x l ) + f{x 2 )- 1 

ri \ _ 2B(A+A)sinh(2Aa) + (B 2 -A 2 -2AA) /r\ 

jy x ) - cnsh 2 f2Ax-i y°) 



Therefore, depending on explicit form of f(x), we get a 
chance to find analytically the zero modes Q n (x) of super- 
charge Q + : 

rt„(x) = exp (x(x))uj n (x) = (6) 

= | cosh(Ax + ) sinh(Ax_)|~ a oj„(x); q + u n (x) = 0. 

The zero modes u> n (x) = r] n (xi)r] n (x2) of operator q + 
are expressed in terms of solutions of one-dimensional 
Schrodingcr equations: 

(7) 



a Slightly different notations for coupling constants are chosen 
here as compared with 1231 



d 1 + /(xi) p?n(xi) = e n r] n (xi). 



Fortunately, for the model under consideration, function 
f{x) is such that this one-dimensional model belongs to 
the class of exactly solvable models [35] , (3D] • It is known 
as Scarf II (hyperbolic Scarf) systems [31], and its spec- 
trum and wave functions were built analytically in terms 
of Jacobi polynomials [3"2"] : 

r) n ( x ) = i"(cosh(2Ax))- A / 2A • (8) 

•cxp[-(S/2A)arctansinh(2Ax)]Pi 7 ' ,9) (isinh(2Aa;)); 

7 = -(A/2X + iB/2X + 1/2); /3 = 7*; e„ = -(A - 2nX) 2 . 
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Therefore, all normalizable zero modes Q n {x) (and their The expression in brackets is transformed by Christoffcl- 

components w„) are known and can be used to find the Darboux formula 22.12.1 [33] : 

wave functions of the Hamiltonian H ^ . 

It follows H], HE] directly from the intertwining rcla- P^f \z 1 )P^^\z 2 ) - P^\zi)P^\z 2 ) = (14) 

tions H]) that the variety of zero modes is closed under _ , \fc ra h. ra _ 1 ^n-i i p(i-P)( \r>{i,P)( \ 

the action of operator H^l In other words, the action of " ^ 2 ~ Zl > k —i ^ m=0 ^ m { - Zl > Fm ^ 2 <>' 

H^ onto Q n gives the linear combination of Q's: , , , , , „ T , . , . , , 

° where the constants h n , k n tor Jacobi polynomials are de- 

n fined in [33] and [32]: 

H^n n {x) = y^c nk ^ k {x), (9) 

trf. h _ Q7+<3+l I>+7+l)I>+/3 + l) . 

fc - U " n — Z (2r l + 7 + / 9 + l)r(™+7+/3 + l)' 

where coefficients c nk form the (TV + 1) x (N + 1) con- h = 2~ n (2 , ra+7+/ ll ! , (15) 

^ " n!(n+7-(-p)! v ' 

stant matrix C. Diagonalization of this matrix, if possible, 



will provide both wave functions and energy eigenvalues Inserting ([13)) and (jl4|) into ([T2[) we obtain: 
of H( 2 \ though not all, in general. 

Calculation of matrix C . From this point and below, *■ x 



(A-n)Pi^\z l )Pi^\z 2 ) 



we'll simplify our presentment by choosing the specific _ •/_ iyi, , i 7 _ ~ \ 
value A = 1/2 in all formulas above. It is more conve- 
nient to look for the explicit form of matrix C in terms of _ 2( » +7)( „. +)3) fc „/ 1 „_ 1 „ n _i j_ ( 7 ,,3) , > ( 7 ,/3) , > 
w„, replacing _ff^ 2 ' by its similarity transform: 2n+-y+/3 fc n _ x 2^m=o h m m l Zl J m l z 2; 



ftW ee cxp {-x)H (2) cxp (+x) = (10) Therefore, the relation (TTT]) takes the form: 

(d 2 + d 2 ) + /(a*) + /(x a ) + 2otanh ^ d+ + 

-2acoth ^ cL - a 2 ; ft(2)w » = ( 2e » ~ a >« ~ 2a (-l)'^o • 

• {A-n)P^\z l )P^\z 2 )- 



2 

JV 



and taking into account, that wqPj (zi)Pj (z 2 ) 
2e„ — a + (11) (— lp'cjj, it can be rewritten as: 



/i( 2 'w„(f) = X) fc=0 c„fcWk(d?) 

Since the functions u n arc factorized onto multipliers 2(n+ 7 )(»+ ) 3) fc„ft„-i v n-i l p (7,/3)/ \ P (y,0)( > 

»?n(»i), ^faa) (see ©), the action of /i^ 2 ) can be sim- 2n+7+,3 fc n _ x Lm=o ^-^ (ZiJ^m \z 2 ) 

plificd: 

/i (2) ^„(x) = 

+ S inhx 1 - S inhx 2 (^coshxi ft - cosh x 2 9 2 J w„. M 2 ) Wn = f 2 e„ - a 2 - 2a(A - roW - (16) 

Introducing new variables zi = is'mhxi; z 2 = ismhx 2 , 4«(n+7) (n+/3) k n h n -i y-m-i (-1)" "*,. 

we consider separately a part ol (|1 lj) : 

Thus, according to definition ([9]) the matrix elements of 
coshxi di - coshx 2 d 2 )oj„ = i(-l) n u • (12) C are: 

Cn, m = 0, /or m > n; c„,„ = 2e„ - a 2 - 2a(A - n)\ 

2afc„h„_i(-l) m j. 



«i - 2 2 j AP^^ijPi 7 ^^) + n(zi,z 2 ) 

The function II appeared in p^|) can be written as: 

Similarly to the models with Morse jS], [H] and Poschl- 

U(zi, z 2 ) = (1 - z 2 )d Zx P^ ,P \zi)P}?'P \z 2 ) - Teller [H] potentials, matrix C is triangular, and its diag- 

_m _ z 2 )p^^h z )Q P^'P\z ) 0na ^ e l ernen t s gi ye immediately the energy eigenvalues of 

H^ : 
and by means of relations 22.8.1 of [33] between Jacobi 

polynomials P n {Zi) and their derivatives over argument, E\ = 2e„— a — 2a(A— n) = — (A— n) — (A— n+a) . (17) 
it can be transformed as follows: __, 

Diagonalization of C . In order to find the correspond- 

II(zi,Z2) = — ra(«i — z 2 )Pn ' (zi)P n ' (z 2 ) + i n g wave functions, one has to diagonalize the matrix C, 

2(n+y)(n+f3) ( D (7,/3)/ n D (7,/3)/ x i- e - to nnd matrices B and A = diag(X a , Ai ,..., Ajv), such 

+ 2 n +7+/3 l^n-1 ^l)^« ( Z 2)- that . 

v JV 

-P^\zi)P^\z2)\ (13) BC = CA; «- ^^c fei = A^. (18) 

/ fc=o 
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The procedure of diagonalization is specific due to triago- 
nality of C, and the corresponding algorithm was proposed 
for the Morse model in |S]. If this task will be solved, 
and matrix B will be found, one can construct a variety 
of eigenfunctions of H'- 2 ' as linear combinations of zero 
modes f2„ : 



,v 



vf 



(2) 



N 



L(£) = E 6 »< ^)- 



(19) 



z=o 



It is convenient to start an algorithm for solution of 
system of linear equations (fl~8| by solving the first line - 
with i = 0. Indeed, taking successively I = TV, TV— 1, ..., I = 
in (|18j). we obtain that boN is arbitrary normalization 
factor, and other &oz are found itcratively: 



i>0(N- 



(JV-l) 



'JON 



bo(N- 



(N-2) 



C N(N- 



JQN 





C N(N 


-1) 






CNN 


-C(JV- 


-1)(JV- 


-1 


) 




C N(N 


-2) 




+ 


CNN 


-C(JV- 


2)(JV- 


2) 




C(jV-l 


)(JV-2 







CNN~ C(jv-l)(JV~l) CNN — C(N_2)(JV-2) 



Analogously, for the second line i = 1 it follows from ([18] 
that &ijv = 0, the element 6i(jv-i) plays the role of arbi 
trary normalization factor, and further 



bi(N-2) = b 



c (W-l)(W-2) 



(N-2) - l(iV-l)c (JV _ 1)cw _ 1) - C(JV _ 2)(w _ 2) ' 

&1HV-3) = 



l(JV-3) 

This procedure can be continued for I = (TV — 3), (TV — 
A),..., I — 0, and after that, for next lines i = 2,3, ...,TV. 
Finally, the obtained matrix B is also triagonal, but with 
vanishing matrix elements lying under the cross diagonal. 
Elements on the cross diagonal are arbitrary normalizing 
factors (boN, &i(iv-i) above). All other elements can be 
compactly written as followqj: 



Om,p Vm,N— m 



r N-p-m , x / 

J2 (r {m) 
i=i 



(20) 



N — m,p 



where the (TV+ 1) triangular matrices labelled t^" 1 \ m 
0, 1, ..., TV arc defined via the matrix elements of C : 



'n,k 



Cn.k 



CN—m,N—m ~ c k.k 



We stress that in ([2U| the expression ( r 



,(m] 



means the 



l-th power of the matrix r^ m \ The repeated index (TV — m) 
is frozen in (|20|) and not summed over. This expression 
allows to write all elements of the m-th line b miP in terms 
of the matrix r^ m ^ and the arbitrary value of the element 
b m ^N-m on the crossed diagonal. These arbitrary values 
can be fixed by the normalization condition for the wave 
functions ^N-n(x) in (H9|. 



e Let us mention the misprint in the upper limit of summation in 
analogous formula (46) in Q5|. 



Shape invariance of the model. The idea of 

shape invariance is one of the most essential new con- 
tributions to the modern Quantum Mechanics made by 
SUSY Quantum Mechanics approachj. It was formulated 
originally [J] in one-dimensional context, and it provided 
the connection between SUSY and exact solvability of the 
model. Namely, all exactly solvable models were shown to 
obey [2] shape invariance, which allow to solve the models 
in a pure algebraic way, without solution of any differen- 
tial equations. This elegant method was generalized jS], 
[15] , [Hi] to two-dimensional case where it leads only to 
partial (quasi-exact) solvability due to many zero modes 
of second order supercharges. 

First shape invariance. It is easy to check that the 
shift of parameter a — > a = a — 1 transforms H^ 2 ' into 

H( 2 \x;h)=H^\x;a), (21) 

where both Hamiltonians ([2]) with potentials ([3]) are in- 
tertwined according to ([1]). As usual |5j, [T5], shape in- 
variance property allows to build the whole tower of new 
wave functions starting from any known wave function 



(22) 



*, ( 1 2 L(x;a)=Q-(a)Q-(a-l)... 
...Q-{a-m + l)^l{x]a-m)] m= 1,2,3,.. 

Herewith, the energy eigenvalues of ^n,m(S; a) are: 

E%l(a) = E^(a-m). 

As a principle wave function & n Q , in (|23|) . one can take an 
arbitrary wave function ^ n (x; a) obtained in the previous 
Section (sec (fH))) with index n) . Its energy E^ is given 
by (fTT|) . and therefore, the energies of states (f2"2)l are: 



(23) 



E?Ua) = -{A -n) 2 -(A-n-m + a) 2 



(24) 



One may verify that for a — —k these energy levels co- 
incide with a part of full spectrum which was found in 

H3J. 

Second shape invariance. It was mentioned in re- 
cent paper [7] that among different solutions of two- 
dimensional SUSY intertwining relations obeying shape 
invariance property, some pairs are equivalent to each 
other up to a linear transformation of coordinates. In par- 
ticular, the model considering in the present paper (the 
model (A5) of [7]) is equivalent to the model (A8) in [7]. 

One can check that the partner Hamiltonians iJW't 2 ) 
on the plane y = (j/i, j/2) with potentials: 



y(ih(2)^ =g 



sinh 2 (Ay2) cosh 2 (Ayi) 



To be honest, it should be noted that this property was al- 
ready known in the framework of well known Factorization Method 
of E.Schrodinger [34] . |29j . but in a slightly different (and not so 
transparent) form. 
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aA(2b=Fl)sinh(Aj/+)-2A 2 b(fc=Fl)+g 2 /2 
cosh 2 (Ay+) 

gA(2b^l) sinh(Ay )-2A 2 fc(b^l)+g 2 /2 _ 
' cosh 2 (Ay_) ' 

y±=yi± y 2 



(25) 



For the particular case of 2 A = 1, the solution of sys- 
tem (|29[) is b = A; a = B, and therefore, the energies in 
([27)) coincide with levels (f24f . By means of laborious but 
elementary calculations, one can check straightforwardly, 
that 

Q~(x;a, b)Q~{x; a, b — 1) = Q~ (x; a,b)Q~ (x; a — 1, 6). (31) 

Similarly to the situation [33] for the Morse model, this 
identity means that starting from the arbitrary principle 
state \&q q (x; a, b) , one may combine both shape invari- 
Q ± {y) — ^d y+ d y _ ± 4C+(y+)dy_ ± AC-{y-)d y+ + B(y); ances (in parameter a, and parameter b) to obtain wave 



are solutions of intertwining relations of the form (JTJ) : 

H^{y)Q+{y) = Q+{y)H^(y) 
with the supercharge^: 



(26) 



B(i7)=4C + (y + )CL(y_H 
C+(y+) = &Atanh(Ay+) 
C-{y-) = 6Atanh(Ay_) 



/i(yi) + /2(y2); 



,(2 



2cosh(Aj/ + ) ' 

a . 

2cosh(At/_) ' 



A(yi) 



cosh- (Ayi) ; 



h(V2) 



ih-(At/ 2 



According to [7], this model is shape- invariant: H^ 2 '(y, b — 
1) = H^ >(y, b), and therefore, if one knows any (principal) 
wave function vD^ ' (y) of if ^ 2 ^ (y) , then the whole tower of 
wave functions can be built due to shape invariancc: 



functions &n,m(x\ a, b). In other words, one can move by 
different paths in the plane (a, b) of parameters. 



E.V.K. is indebted to the non-profit foundation "Dy- 
nasty" for financial support. 
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*$»(# b) = Q-(y; b)...Q-(y; b-m + l)^ (y; b - m) 
m = 1, 2...; E ni . m (b) = E nfi (b - m). 

At the first sight, the model ([23]) has no relation to the 
Scarf II model ([3]) . But as it was noticed in [7] , they are 
related by substitution of coordinates: x+ = j/i; X- =2/2, 
since (l25l) can be rewritten as: 



(28) 



mW(Z)=g(- 



^ sinh 2 (Aa;_) cosh 2 (Aa:+) ' 

aA(2b^l)sinh(2Az 1 )~2A 2 b(b^l)+a 2 /2 . 
cosh 2 (2Aa;i) ' 

. ctA(2b=Fl) sinh(2Aa;2)-2A 2 b(b^l)+a 2 /2 
' cosh 2 (2Aa;2) ' 



If one identify the coupling constants in 
follows: 



[9] 



and © as 
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-- A 2 a(a + 1) 
* 2 -4A 2 6(fo 



aA(26+l) = B{A + \y 7 (29) 
1) = (B 2 -A 2 - 2AX), 



it becomes evident, that the second supcrpartners are pro- 
portional: H^ 2 \x) = 2H i - 2 ^(x), and thus, the Hamiltonian 
H^> (x) satisfies simultaneously two intertwining relations 
- with H (1 ~> and #« : 

iJ (1) (.f; a, 6)-Kff (2) (.f; a, b) = 2H {2 \x; a, b)~2H (1) (.f; a, b), 

(30) 
where the sign -=- means intertwining, and both shape in- 
variance parameters a, b of potentials are written explic- 
itly 



J We use everywhere in the present paper the same definitions of 
coefficient functions C± as in [7], [23] ■ They differ from definitions 
of C± in all earlier papers by the factor 1/4. 
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